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The quantum interference and orbital filling effects on the thermoelectric (TE) properties of 
quantum dot molecules with high figure of merit are illustrated via the full solution to the Hubbard- 
Anderson model in the Coulomb blockade regime. It is found that under certain condition in the 
triangular QD molecule (TQDM), destructive quantum interference (QI) can occur, which leads 
to vanishing small electrical conductance, while the Seebeck coefficient is modified dramatically. 
When TQDM is in the charge localization state due to QI, the Seebeck coefficient is seriously 
suppressed at low temperature, but highly enhanced at high temperature. Meanwhile, the behavior 
of Lorenz number reveals that it is easier to block charge transport via destructive QI than the 
electron heat transport at high temperatures. The maximum power factor (PF) in TQDM occurs 
at full-filling condition. Nevertheless, low-filling condition is preferred for getting maximum PF in 
serially coupled triple QDs in general. In double QDs, the maximum PF can be achieved either 
with orbital-depletion or orbital-filling as a result of electron-hole symmetry. Our theoretical work 
provides a useful guideline for advancing the nanoscale TE technology. 


I. INTRODUCTION 

To design solid-state coolers and power generators,[1-7] 
many efforts seek efficient thermoelectric (TE) materials 
with the figure of merit ( ZT ) larger than 3, which will 
lead to efficient TE conversion, making the TE device 
competitive with conventional air conditioners and power 
generators. [1,2] Furthermore, the energy harvesting by 
using TE materials is one of several advanced techniques 
considered for clean energies. The optimization of ZT = 
S 2 G e T/n depends on the electrical conductance (G e ), 
Seebeck coefficient (S), and thermal conductance ( k ). T 
is the equilibrium temperature. These physical quantities 
are usually related to one another. Mechanisms leading 
to the enhancement of power factor (PF = S 2 G e ) would 
also enhance the thermal conductance. Consequently, it 
is difficult to obtain ZT above one in conventional bulk 
materials. [1] 

Recently, quantum dot superlattice (QDSL) nanowires 
with impressive ZT values (larger than one) have been 
demonstrated experimentally. [8] The power factor and 
thermal conductance become independent thermoelec¬ 
tric variables under the condition K e /n p h <C 1, where 
K e and K p h denote, respectively, the electron thermal 
conductance and phonon thermal conductance. [8] In the 
Coulomb blockade regime, electron transport process is 
seriously suppressed by the electron Coulomb interac¬ 
tions, causing K e and G e to be reduced simultaneously. [9] 
On the other hand, under the condition K e /n p h "C 1, 
one can increase the power factor (PF) [10-15] and de¬ 
crease the phonon thermal conductance [16-21] simulta¬ 
neously to optimize ZT. The ZT enhancement of QDSL 
mostly arises from changes induced by dimensional con¬ 
finement in the electronic band structures as well as en¬ 


hanced phonon scattering resulting from the scattering 
at nanowire surface and interfaces surrounding QDs. [1,2] 


Although many theoretical efforts have investigated 
the PF enhancement of QDs[10-15] and molecule 
junctions[22-24], quantum interference (QI) and orbital 
filling effects on the PF optimization are still puzzling 
due to the difficulty to treat many body effect reliably 
in either QDSLs or molecules. [25] For example, in the 
triangular QD molecule (TQDM), there are 923 electron 
correlation functions appeared in the 4752 Green’s func¬ 
tions need to be solved, when electron Coulomb interac¬ 
tions are turned on, even though one just considers one 
energy level for each QD. As a consequence, a theoretical 
framework to treat adequately the many-body problem of 
a molecular junction remains elusive. QI is a remarkable 
effect, which influences the charge transport of QD array 
[6-8] and molecules [22-24] with multiple quantum paths. 
When the energy levels of QDSLs or molecules are below 
the Fermi energy of electrodes, the orbital filling effects 
can not be avoided for charge/heat transport. To reveal 
the many-body effects on the TE properties of QD junc¬ 
tion system, we theoretically investigate QD molecules 
(QDMs), including double QD molecule (DQD), serially 
coupled triple QDs (SCTQD), and TQDM with a full 
many-body solution which takes into account all correla¬ 
tion functions resulting from electron Coulomb interac¬ 
tions in QDMs. Our results indicate that the effects of 
QI and structure-dependent orbital filling are quite sig¬ 
nificant in determining the behaviors of TE coefficients 
for junctions of QDSL or molecules ( such as benzene). 
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II. FORMALISM 

Here we consider nanoscale semiconductor QDs, in 
which the energy level separations are much larger than 
their on-site Coulomb interactions and thermal energies. 
Thus, only one energy level for each quantum dot needs 
to be considered. An extended Hubbard-Anderson model 
is employed to simulate a QDMs connected to electrodes 
(see the inset of Fig. 1(b)). The Hamiltonian of the 
QDM junction is given by H = H 0 + Hqjj: 

Ho = £fcQfc. cr Qfc,a-+fc.o-fofc.a (i) 

k,cr k,cr 

^ ^ VkiLd^ pCLkiCr H” ^ ^ C.C 

k,cr k,cr 

where the first two terms describe the free electron gas of 
left and right electrodes. a\. a (b j) CT ) creates an electron of 
momentum k and spin a with energy e in the left (right) 
electrode. 14y (£ = L, R) describes the coupling between 
the electrodes and the left (right) QD. d\ a (c4 jfT ) creates 
(destroys) an electron in the £-th dot. 


H QD = £ E e n ejCT + ^2 Uini ta ni,g (2) 

£,cr t 
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where Et is the spin-independent QD energy level, 
and ri( ttT = d\ a de,, a . Notations Ue and Utj describe 
the intradot and interdot Coulomb interactions, respec¬ 
tively. t( j describes the electron interdot hopping. Not¬ 
ing that the interdot Coulomb interactions as well as in¬ 
tradot Coulomb interactions play a significant role on the 
charge transport for semiconductor QDs [9] and molecu¬ 
lar structures. [22-24] 

Using the Keldysh-Green’s function technique, [26,27] 
the charge and heat currents from reservoir a to the 
QDM junction are calculated according to the Meir- 
Wingreen formula 

J* = J E/ deTJ(e)[G< a (e) + f a (e)(G r j<7 (e) 

jo 

- G a oa m (3) 

Qa = de(e-f, a )r«(e)[G< (e)/ a (e) 

jo- 

(G r ja (e)-G<] a (e))], (4) 

Notation T“ = |14,a,f | 2 ^(e — eQ is the tunneling rate 

between the a-th reservoir and the i -th QD. f a (e) = 
l/{exp[(e — Ha)/kBT a ] + 1} denotes the Fermi distribu¬ 
tion function for the a-th electrode, where p, a and T a are 
the chemical potential and the temperature of the a elec¬ 
trode. fiL — M R = AU and 74 — Tr = A T. e, h, and &b 


denote the electron charge, the Planck’s constant, and 
the Boltzmann constant, respectively. G^ a (e), Gj a (e), 
and Gj a (e) are the frequency domain representations of 
the one-particle lessor, retarded, and advanced Green’s 
functions Gf a {t,tf) = i(d) a (t')d h(J (t)), G r ja (t,t') = 
-iO(t - <7 (*')}). and G a i<T {t,t') = iQ(f - 

t)({dj,<r{t)>dj (t')}), respectively. These one-particle 
Green’s functions are related recursively to other Green’s 
functions and density-density correlators via the few- 
body equation of motion, [25] which we solve via an itera¬ 
tive numerical procedure to obtain all n-particle Green’s 
functions and correlators for the QDM. The coupling be¬ 
tween electrons in the leads and electrons in QDM is 
included in the self-energy term, T“. Because of this ap¬ 
proximation, our procedure is valid only in the Coulomb 
blockade regime, but not the Hondo regime. [28] 

Thermoelectric coefficients in the linear response 
regime are 
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The quantity sf ( t ) ' obtained by solving the varia¬ 
tion of the equation of motion with respect to the change 
in Fermi-Dirac distribution, f a (e). Here we have assumed 
the variation of the correlation functions with respect 
to Sf a (e) is of the second order. Note that we have to 
take AU — > 0 for the calculations of (j^)at=o an d 
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(j 2 $v)at=o- Meanwhile we taken AT —> 0 for Eqs. (9) 
and ( 10 ). Ep is the Fermi energy of electrodes. There is 
a Joule heating term AV x J a arising from Eq. (4) which 
can be ignored in the linear response regime. 

It is a challenge to prove that whether the Eqs. (3) and 
(4) can be expressed in terms of transmission coefficient 
(or Landauer’s expression) when one takes into account 
all Green’s functions and correlation functions fully. To 
gain deeper insight into the electron correlation effect, we 
adopt the approximation procedure of Ref. [9] and obtain 
G e = e 2 C 0 , S = -Ci/(eTC 0 ) and K e = P(C 2 - C\/C Q ), 
which can be calculated by a closed form expression for 
the transmission coefficient. C n is given by 

Cn =U deT LR (e)(e - EP)^, ( 12 ) 

where Tlh(c) is the transmission coefficient, /(e) = 
l/(exp( e ~ EF ^ kBT + 1). The thermoelectric coefficients 
determined by Eq. ( 12 ) are employed to compare with 
results of Eqs. (5)-(7). Because the condition K e /n P h -C 1 
is readily satisfied in the Coulomb blockade regime,[9] the 
optimization of ZT can be improved by finding the best 
power factor (PF = S 2 G e ). 

III. RESULTS AND DISCUSSION 
A. Effect of destructive QI 

Firstly, we consider TQDM with electron hopping 
strengths tpc = tcR = f c / thR- Figure 1 shows the 
electrical conductance and Seebeck coefficient as func¬ 
tions of the energy level of central QD (Ac = Ec — Ep) 
for various coupling strengths, tLR at k B T = ir 0 . In 
the small Ac regime (Ac < 10r o ) G e is found to be 
insensitive to tLR, because the charge transport is dom¬ 
inated by the near-resonant tunneling process through 
levels at Ep + y/2t c , Ep and Ep — \[2t c . On the other 
hand, G e changes significantly with the variation of tLR 
in the regime of Ac/r 0 > 10. For tp R = O.ir 0 , G e is 
vanishingly small at Ac 4 = 90Fo- The minimum of G e 
occurs at lower Ac if t.pR is increased. For example, the 
minimum of G e occurs at Ac 3 = 45Fo, Ac 2 = 30Fo, 
and Aci = 22.5r 0 for tLR = 0.2r o , t B R = 0.3r o and 
tLR = 0.4r 0 , respectively. This behavior is caused 
by the quantum interference (QI) between two paths. 
This direct coupling (t lr ) between the outer QDs pro¬ 
vides one path. The long distance coherent tunneling 
(LDCT)[25] mediated by the central QD provides an¬ 
other path, which leads to an effective electron hopping 
strength t e ff = —tLctcR /between the outer QDs. 
Such an effective coupling between outer dots can be used 
to manipulate the spin entanglement between dots sepa¬ 
rated by a long distance. [29-31] We note that the vanish¬ 
ingly small G e occurs when the condition of \t e ff\ = tLR 
is met. To depict the behavior of Fig. 1(a), we also cal¬ 
culate G e by using the approximate expression, 7/ fl (e) 


for transmission coefficient as given below. 

TUe) (13) 

_ L C rPiFq! _ 

|MlM2/i3 - t 2 CR p 1 - t 2 LC p 3 - t 2 LR p, 2 - 2t L RtLCtcR\ 2 ’ 

where = e - E L + iT L , = e — Ec and p 3 = 
e — Er + ir#. Pi denotes the probability weight of elec¬ 
tron transport through TQDM in an empty state, which 
is determined by the one particle occupation number and 
on-site two particle correlation functions resulting from 
electron Coulomb interactions. Pi equals to one in the 
absence of electron Coulomb interactions. The numera¬ 
tor Fqi = (t lc^cr + tLRH 2 ) 2 causes the QI effect. The 
first term and the second term of Fqi describe, respec¬ 
tively, the upper path through the central QD and lower 
path for direct hopping between the outer QDs. Due to 
t lc = tcR = / / tLR , there are three poles in the de¬ 
nominator of Eq. (13). If we take F = 0, the poles occur 
at 

e± = E ° + E ° + tLR ± \^{E C -Eo- t LR ) 2 + 8 t 2 c 

£0 = Eq — tLR- (14) 

Here, Eq = El = Er = Ep in Eq. (14). For higher 
symmetry TQDM (t LR = t L c = tcR = t c and E e = E 0 ), 
we have e+ = Eq + 2 t c , e_ = Eq — t c and eo = Eq — t c . 

Once Ac = Ec — Ef 2t c , the lowest energy level 
is given by eo = Eq — tLR- Keeping only the resonant 
channel at eo, we obtain 

r< - 2e2 r?r ( t e ff + t L R) 2 _ Pi 

e ~ h Ak B T (t LR + t eff ) 2 + (r/2) 2 cosh 2 ($£,)’ 

(!5) 

where t e ff = — t 2 /Ac- This expression can well explain 
the destructive QI behavior shown in Fig. 1(a). Accord¬ 
ing to Eq. (13), the QI observed in Fig. 1(a) will dis¬ 
appear as tLR = 0 (the case of SCTQD). In Sanchez et 
al[32] it was proposed that destructive QI of two superex¬ 
change trajectories can also occur in SCTQD. Note that 
such a effect is easily masked by the background current 
at finite temperatures. Thus, it is not as robust as the 
QI effect described here. Using Eq. (12) and considering 
the eo pole of Eq. (14), we can show that S = tLR/T is 
positive and independent on Ac- However, S = tL R /T 
can not describe the behavior of S shown in Fig. 1(b) 
(obtained by the full calculation), which is a function of 
Ac. This implies that the resonant channels involving 
Ec can not be ignored for calculating S in the regime 
of t e ff /tLR > I- The positive sign of S indicates that 
the hole diffusion dominates over electron diffusion at 
low temperature (k B T = lTo). Here holes are defined 
as missing electrons in states below Ef. The results of 
Fig. 1(b) show that the Seebeck coefficient is seriously 
suppressed under the destructive QI effect. Such a behav¬ 
ior is quite different from the general behavior of thermal 
voltage observed in QD junction systems. [10-15] (See re¬ 
sults of Figs. 5-8 and note that the minimum of S occurs 
at the maximum of G e ). 
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To further clarify the above destructive QI behavior, 
we plot the electrical conductance and Seebeck coeffi¬ 
cient at various temperatures for t^R = 0.3To in Fig. 2. 
It shows that G e is suppressed with increasing tempera¬ 
ture. Meanwhile, the destructive QI effect on G e is very 
robust with respect to temperature variation when kpT 
is tuned up to 5Fo. The results of Fig. 2(a) imply that 
a single-electron QI transistor is achievable even at room 
temperature. [33] The change of sign for S from positive 
to negative indicates that the resonant channels above 
Ep become important with increasing temperature with 
electron contribution dominating over hole contribution. 
In particular, we noticed a large enhancement of the neg¬ 
ative S peak at Aqj (the value of Ac where QI occurs) 
as temperature increases. Namely, the thermal voltage 
changes sign and enhances with increasing temperature. 
We also noticed that peak value of |S| reaches a maxi¬ 
mum near kpT = 4To and decreases afterwards. 

To reveal the importance of electron correlation effects 
and understand the interesting behavior of S in Fig. 2(b), 
we recalculate the results of Fig. 2 with the procedure of 
Ref. [9]. There are 32 configurations in the transmission 
coefficient of Eq. (12) for electrons with spin a in the 
electrodes. The expression of Eq. (13) is for an empty 
TQDM. The curves in Fig. 3 have one-to-one correspon¬ 
dence to those of Fig. 2. We found significant differences 
between Fig. 2(a) and Fig. 3(a), when Ac is smaller than 
20 To. This implies that the electron correlation effects 
become important when Ec closes to Ep. In particu¬ 
lar, two-particle interdot correlation functions. However, 
both approaches give similar QI effect with the same 
value for A qj. This is because QI is mainly caused by 
the single-particle transport process. Fig. 3(b) also shows 
enhancement of |Sj at A qj with increasing temperature, 
but the degree of enhancement is overestimated com¬ 
pared with the full calculation. The temperature depen¬ 
dence of S max (at Aqj) is plotted in the inset of Fig. 3(b). 
It is seen that S max obeys the simple relation S ~ 
when kpT > 2To. One can show that the next impor¬ 
tant contribution to S is from the resonant channel near 
e = Eq + Ulr — t pr, but not from e + given by Eq. (14). 
Note that the contribution of this resonant channel to G e 
is small compared to the eo channel. To demonstrate the 
effect of electron Coulomb interactions, we plot the curve 
with red triangle marks for Ug = Ug.j = 0 at = 2To. 
Now, Smax = tLFt/T becomes very small in the absence 
of electron Coulomb interactions. This proves that the 
enhancement of \S m ax\ results from the resonant chan¬ 
nels involving electron Coulomb interactions. This also 
explains the enhancement of |iS maa; | shown in Fig. 2(b). 

According to the results of Figs. (1) and (2), the de¬ 
structive QI effect can blockade the charge transport. 
Next, we clarify how QI influences electron heat trans¬ 
port. Fig. 4 shows the electron thermal conductance (n e ) 
and Lorenz number ( L z = n e /(G e T)) as functions of 
Ac. The physical parameters adopted in the calculation 
are the same as those for Fig. 2. Like G e , the electron 
heat transport can be modulated by tuning the central 


QD energy level. We see that the electron thermal con¬ 
ductance is also vanishingly small at Aqj = 30Fo due 
to the destructive QI. There is the Wiedemann-Franz 
law in metals to link two physical quantities ( G e and 
K e ). Here, we recheck the Wiedemann-Franz law by 
calculating the Lorenz number in Fig. 4(b). A highly 
temperature-dependent behavior of L z occurs at Aqj. 
The fact that L z larger than one indicates that it is eas¬ 
ier to block charge transport via destructive QI than the 
electron heat transport. In metals the Lorenz number is a 
universal constant L z = ^-(/cg/e) 2 , which is independent 
of T. Our results obviously violet the Wiedemann-Franz 
law near A c = Aqj for all temperatures. The viola¬ 
tion of Wiedemann-Franz law is a typical phenomenon 
for nanostructures with discrete energy levels. [3] The 
QI effects on the TE properties are illustrated by us¬ 
ing QDMs with N=3 as an example. The asymmetrical 
features of tLC = tcR ^ t lr , Ulc = Ucr ^ Ulr and 
El = Er ^ Eq considered in TQDM can be used to re¬ 
veal the effects of of QD size fluctuation and nonuniform 
interdot separation on the charge transport of QDSL sys¬ 
tem. [6,8] 


B. Orbital filling effects 

Many theoretical studies have used the first-principles 
method (within mean-field approximation) to predict the 
TE coefficient of realistic molecules with orbital below 
the Ep of electrodes. [22-25] Such a mean-field approach 
does not take full account of the electron correlation ef¬ 
fects. Therefore, the predicted G e , S and PF can not 
reveal the realistic TE properties of molecule junction 
systems in the Coulomb blockade regime. To reveal the 
orbital charge-filling (orbital below Ep) effects on the 
thermoelectric properties of TQDM, we plot the total oc¬ 
cupation number ( N = ^^Nl^+Nc^+Nr^)), G e and 
S as functions of gate voltage V g (Eg = Ep + 30r o — eV g ) 
at T = 0.3r o , tgj = 3r 0 and Ugj = 30r o for various 
temperatures in Fig. 5. Here, the QD energy levels are 
placed at 30r o above Ep when eV g = 0. The average to¬ 
tal occupation number (N) display an staircase behavior 
with six plateaus. The step edges are broadened as tem¬ 
perature increases. The height of each staircase is equal 
to one. The onsets of these plateaus are roughly deter¬ 
mined by Ugj and Ug. Therefore, the electron number of 
TQDM is tuned from one to six with increasing gate volt¬ 
age. The spectrum of G e exhibits a Coulomb oscillation 
behavior with respect to gate voltage. Six main peaks 
of G e labeled by e n correspond to the particle-addition 
energy of TQDM with different many-body states. For 
kpT = lTo, there are several insulating states (with van¬ 
ishingly small G e ). Such insulating states are not due 
to the QI effect, but the Coulomb blockade effect. The 
regime of vanishingly small G e is related to a very small 
transmission coefficient resulting from the absence of res¬ 
onant levels near Ep. The insulating state of a single 
benzene molecule calculated by the DFT method was 
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suggested from the destructive QI effect. [22-24] The sug¬ 
gestions from Refs. [22-24] are unclear due to the de¬ 
ficiency of the mean-field theory (DFT). Although the 
behavior of G e with respect to kpT is a typical behavior 
that the magnitude of G e peak is suppressed with in¬ 
creasing temperature, meanwhile the width of each peak 
becomes broadened, we find a more dramatic reduction 
of eg peak when kpT is tuned from lTo to 2To. In ad¬ 
dition, the peak position is shifted with increasing tem¬ 
perature. Several secondary peaks resulting from excite 
states around these main peaks are washed out by in¬ 
creasing temperature. The separation between £3 and £4 
is the so called ” Coulomb gap” (these two peaks are sepa¬ 
rated by the intradot Coulomb interaction Ue).[ 25] Using 
the middle Coulomb gap as a reference point, the G e 
spectrum does not maintain the mirror symmetry. This 
is different from that of DQDs. 

The Seebeck coefficient shown in Fig. 5(c) also shows 
an oscillatory behaviors with respect to gate voltage. 
When QD energy levels are above Ep the Seebeck coef¬ 
ficients are negative, which indicates that electrons dom¬ 
inate the diffusion process of thermoelectric properties. 
When QD energy levels near Ep, S almost vanishes due 
to the electron-hole balance. The positive S indicates 
that holes through resonant level below Ep become ma¬ 
jority carriers. Such bipolar effects are the interplay 
between holes and electrons. When the resonant lev¬ 
els corresponding to the two-electron and three-electron 
states of TQDM (£2 and £ 3 ) are near Ep, once again the 
Seebeck coefficients become very small. We note that 
when G e reaches the maximum, the Seebeck coefficient 
(S = AV/AT) reaches zero. This is different from that of 
Fig. 2. The Seebeck coefficient near the onsets of plateaus 
for N = 1, 3,4, and 6 are highly enhanced with increasing 
temperature. Like the G e spectrum, the Seebeck coeffi¬ 
cient does not exhibit a mirror symmetry. We find promi¬ 
nent secondary oscillatory structures of S (clearly notice¬ 
able at kpT = lTo), which arise from the secondary res¬ 
onant levels around the main peaks of G e . Compared 
to the spectrum of G e at low temperature, the spectrum 
of S can clearly resolve the structures arising from the 
excited states (secondary resonances) of TQDM. Conse¬ 
quently, the measurement of S' is a powerful means to 
reveal various configurations of TQDM arising from the 
many-body effect. Unlike the case of metallic QDs,[34] 
the spectrum of S does not show the periodically oscillat¬ 
ing oscillatory structure with respect to the gate voltage 
due to U(; ^ U(j and electron correlation effects. The in¬ 
terdot Coulomb interactions as well as intradot Coulomb 
interactions play significant roles for electron transport. 

To examine the interdot Coulomb interaction effects, 
Figure 6 shows (a) the electrical conductance, (b) See¬ 
beck coefficient and (c) power factor as functions of V g 
for three different values of Ut,j at kpT = lTo- Other 
physical parameters are the same as those of Fig. 5. From 
the results of G e , we find that the interdot Coulomb 
interactions not only influence the peak positions, but 
also change the magnitude of each peak. For Ufj = 0, 


the G e spectrum shows only four peaks (see blue dashed 
line). The first two peaks (separated by 3t c ) correspond 
to the resonant levels £bd = Eq — t c (bonding state) and 
cab = E 0 + 2t c (antibonding state). It is noted that the 
first G e peak at epD is significantly reduced and the sec¬ 
ond peak at cab is totally suppressed when Uej becomes 
finite. 

The peak structures of the G e spectrum for Ufj = 
15Fo is quite similar to that for U( t] = 30Fo, except 
the spacings between peaks are different, which leads to 
very different behavior in the Seebeck coefficient. Pre¬ 
vious theoretical works predicted the thermal power (or 
Seebeck coefficient) spectra without considering interdot 
Coulomb interactions. [15]. This is inadequate for study¬ 
ing semiconductor QD molecules with nanoscale sepa¬ 
ration between QDs. Our calculations indicate that to 
achieve large maximum power factor (PF ma:r ) in TQDM 
with homogenous QD energy levels and electron hop¬ 
ping strengths, the full-filling condition (with six-electron 
state) is preferred. This is due to the fact that three 
resonant channels associated states localized at dot L, 
C, and R (with energies at Ep + Uq + 2 Ulc + 2 Ulr, 
Ec + Uq + 2 Ulc + 2 Uqr, and Ep + Uq + 2Ulr + 2Ucr, 
respectively) are all aligned. It is interesting to see what 
would happen in other QDM structures when the reso¬ 
nant conditions are changed. 

Figure 7 shows the electrical conductance (G e ), See¬ 
beck coefficient (S) and power factor (PF) of DQD 
as functions of V g at different temperatures. Physi¬ 
cal parameters are the same as those for Fig. 5. For 
DQD, the four G e peaks corresp ond to = Eq — tpR , 
£2 = Eq + tLR + 1 sATLlr — i^/({7o — Ulr) 2 + 

£3 = Eq - t LR + UsAEIlr + ^(Uo-UlrY + ltii R 
and £4 = Eq + Uq + 2Ulr + tLR■ The G e spectrum 
shows a mirror-symmetry behavior with respect to the 
middle of the Coulomb gap. Due to the symmetry of 
DQD structure, the PF spectrum also shows the mirror 
symmetry. The maximum PF occurs at either orbital- 
depletion ( TV* < 1) or orbital-filling (N = 4) condition for 
DQD. Note that unlike G e and PF, S does not show the 
mirror symmetry. 

Figure 8 shows G e , S and PF of SCTQD as functions 
of V g at different T. For simplicity, we adopt tLR = 0 and 
Ulr = 0 , since they are typically much smaller than the 
corresponding parameters related to the central dot. The 
G e spectrum in Fig. 8 does not show the mirror symme¬ 
try as a result of the inhomogeneous interdot Coulomb 
interactions ( Ulc — Ucr ^ Ulr)- The peaks labeled £jv 
(N = 1 , • • • ,6) correspond to the transitions where the 
total particle numbers in the SCTQD changes from IV — 1 
to N with £4, £2, and £4 being more prominent. The ei 
peak comes from electrons tunneling through the empty 
SCTQD via the resonant level e± = E 0 — \[ 2 t c . The £2 
peak ( £ 2 = EQ+\/ 2 t c +U lr) corresponds to charge trans¬ 
port in the presence of another electron localized in one 
of the two outer dots. (Here we have chosen Ulr = 0) 
This is verified by examining Nq^, which is negligible 
until eV g is near £3. £4 = Eq + Uq + Ulc{Ucr) corre- 
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sponds to the injection of an electron into the left (right) 
dot of SCTQD with all three dots each filled with one 
electrons, while €5 = Eg + Uq + Ulc(Ucr) + 2t e ff corre¬ 
sponds to the addition of another electron in either left or 
right dot to the configuration associated with £4, where 
teff = t L ctcil/ {Ef — (Ec + 2 Ulc + 2 Ucr )) = 0.9r o 
is the effective hopping strength mediated by the center 
dot. Therefore, the splitting between £4 and £5 can be 
observed by increasing tpj ■ The Seebeck coefficient spec¬ 
trum can reveal fine structures not resolved in the G e 
spectrum. For example, the £3 and £6 peaks are hardly 
noticeable in G e but clearly observable in S. The peak 
at £3 = E 0 + Urc + Ucr is caused by the charging of 
the center dot with two outer dots each filled with one 
electron and £6 describes the charging of SCTQD with 
a second electron into the center QD with two outer 
dots each filled with two electrons. In strong contrast 
to TQDM, both G e and PF of SCTQD are very small at 
the maximum filling condition (IV = 6), which occurs at 
eV g = 250r o where Ec + Uo + 2Ulc + 2 Ucr is aligned 
with Ep. This is because the three resonant channels as¬ 
sociated states localized at dot L, C, and R (with energies 
at El + Uq + ZUlc + ^Ulr, Ec + Uq + 2Ulc + 2Ucr , and 
Er + Uo + 2Ulr + 2Ucr, respectively) are misaligned, 
since Ulr = 0 (or Ulr -C Ucr in general). Conse¬ 
quently, charge transport is seriously suppressed for QD 
molecule with orbital-filling condition, which is expected 
to be a general feature in serially coupled QDMs with 
more dots. For the present SCTQD, the maximum PF 
occurs at the combined £4 and £5 peak as shown in Fig. 8. 
In this case, the charge transport is enabled via the long 
distance coherent tunneling mechanism, which means 
that the localized electron of the left dot can be trans¬ 
ferred to the right dot via an effective hopping strength 
mediated by the central dot.[25] When we consider Ulr 
or tune tgj away from 3Fo (not shown here), the maxi¬ 
mum PF for SCTQD would occur at the orbital-depletion 
condition. 


IV. SUMMARY 


system, we have theoretically investigated the TE coef¬ 
ficients of TQDM with a full many-body solution. Our 
theoretical work can serve as a guideline for the design of 
nanoscale TE devices. The destructive QI of G e is found 
to be very robust with respect to temperature variation. 
This implies that it is possible to achieve the destruc¬ 
tive QI of G e at high temperatures. Furthermore, the 
Seebeck coefficient can be made large with a vanishingly 
small electrical conductance when temperature increases. 
From the results of various QDMs with a few electrons, 
we have demonstrated that the measurement of Seebeck 
coefficient, S' is a powerful tool to reveal the many-body 
effect of nanostructures, since S is influenced much more 
than G e .To achieve maximum power factor in the case 
of TQDM the orbital-filling situation is preferred (with 
N= 6). However, for SCTQD in general it is preferable 
to the orbital-depletion situation with (IV < 1) and likely 
for serially coupled QDMs with more dots. It is desirable 
to compare current results including all Green functions 
and electron correlations with those calculated by the 
Hartree-Fock approximation (HFA) [22-24] to examine 
the validity of mean-field approach. Such a comparison 
for SCTQD has been carried out. [35] We found that HFA 
works well only for the orbital-depletion situation and 
high temperatures. 
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Ac /r O 


FIG. 1: (a) Electrical conductance (G e ) and (b) Seebeck co¬ 
efficient (S) of TQDM as a function of central QD energy 
(Ac = Ec — Ef ) for different tLR strengths at El = Er = 
Ef, t lc = tcR = t c = 3r 0 and ksT = lFo. We assume 
Ulc = Ucr = 30r 0 , Ulr = 10r o , Ui = Uo = lOOro and 
t l = r iJ = r = o.3r 0 . 



FIG. 3: The curves of Fig. 3 are one to one corresponding to 
those of Fig. 2. The calculation of Fig. 3 only considers the 
single particle occupation numbers and on-site two particle 
correlation functions. 



AC /r 0 


FIG. 2: (a) Electrical conductance (G e ) and (b) Seebeck co¬ 
efficient (S) of TQDM as a function of central QD energy 
(Ac = Ec — Ef) at tLR = 0.3Fo for different temperatures. 
Other physical parameters are the same as those of Fig. 1. 



Ac /r o 


FIG. 4: (a) Electron thermal conductance (ft e ), and (b) 

Lorenz number(K e /(GeT)) as a function of central QD energy 
for different temperatures. The curves of Fig. 4 correspond 
to those of Fig. 2. fco = k\/e 2 . 
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FIG. 5: (a) Total occupation number, (b) electrical conduc¬ 
tance (G e ), and (c) Seebeck coefficient (S) as a function of QD 
energy (Ee = Ef + 30Fo — eV g ) at te,j = 3Fo and Ue,j = 30Fo 
for different temperatures. Other physical parameters are the 
same as those of Fig. 1. 



FIG. 6: (a) Electrical conductance (G e ), (b) Seebeck coef¬ 
ficient (S) and (c) power factor as a function of QD energy 
(Ee = Ef + 30Fo — eV g ) at te,j = t c = 3Fo and ksT = lFo 
for different interdot Coulomb interactions. Other physical 
parameters are the same as those of Fig. 5. 



FIG. 7: Electrical conductance (G e ), Seebeck coefficient (S) 
and power factor of DQD as functions of QD energy (Ee = 
Eq = Ef + 30r 0 — eV g ) for different temperatures. Other 
physical parameters are the same as those of Fig. 5. 



eV g /r 0 


FIG. 8: Electrical conductance (G e ), Seebeck coefficient (S) 
and power factor of SCTQD as functions of QD energy (Ee = 
Ef + 30Fo — eV g ) for different temperatures. Ulr = 0 and 
t^R = 0. Other physical parameters are the same as those of 
Fig. 5. 






































